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Abstract. We report the results of an observation of Cygnus X-1 performed on June 7-11 2003 with INTEGRAL that we 
combine with simultaneous radio observations with the Ryle telescope. Both spectral and variability properties of the source 
indicate that Cygnus X-1 was in an Intermediate State. The INTEGRAL spectrum shows a high-energy cut-off or break around 
100 keV. The shape of this cut-off differs from pure thermal Comptonisation, suggesting the presence of a non-thermal compo- 
nent at higher energies. The average broad band spectrum is well represented by hybrid thermal/non-thermal Comptonisation 
models. However, models with mono-energetic injection, or models with an additional soft component are favoured over stan- 
dard power-law acceleration models. During the 4 day long observation the broad band (3-200 keV) luminosity varied by up 
to a factor of 2.6 and the source showed an important spectral variability. A principal component analysis demonstrates that 
most of this variability occurs through 2 independent modes. The first mode consists in changes in the overall luminosity on 
time scale of hours with almost constant spectra (responsible for 68 % of the variance) that are strikingly uncorrelated with the 
variable radio flux. We interpret this variability mode as variations of the dissipation rate in the corona, possibly associated with 
magnetic flares. The second variability mode consists in a pivoting of the spectrum around ~10 keV (27 % of the variance). 
It acts on a longer time-scale: initially soft, the spectrum hardens in the first part of the observation and then softens again. 
This pivoting pattern is strongly correlated with the radio (15 GHz) emission: radio fluxes are stronger when the INTEGRAL 
spectrum is harder We propose that the pivoting mode represents a 'mini' state transition from a nearly High Soft State to a 
nearly Low Hard State , and back. This mini-transition would be caused by changes in the soft cooling photons flux in the 
hot Comptonising plasma associated with an increase of the temperature of the accretion disc. The jet power then appears to 
be anti-correlated with the disc luminosity and unrelated to the coronal power. This is in sharp contrast with previous results 
obtained for the Low Hard State, suggesting a different mode of coupling between the jet, the cold disc, and the corona in 
Intermediate States. From this interpretation we also infer that the bolometric luminosity jumps by a factor of about 2 during 
the transition hard to soft, suggesting a radiatively inefficient accretion flow in the Low Hard State. 

Key words. Gamma-rays: observations -Black hole physics - Radiation mechanisms: non-thermal -X-rays: binaries; radio 
continuum: stars - X-rays: individual: Cygnus X-1 



1. Introduction 

Cygnus X-1 is the prototype of black hole candidates. Since its 
discovery in 1964 (Bowyer et al. 1965), it has been intensively 
observed by all the high-energy instruments, from soft X-rays 
to y-rays. It is a persistent source most often observed in the 
so-called Low Hard State (hereafter LHS), characterised by a 
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relatively low flux in the soft X-rays (~ 1 keV) and a high flux 
in the hard X-rays (~ 100 keV). In the LHS, the high-energy 
spectrum can be roughly described by a power-law with spec- 
tral index F varying in the range 1.4-2.2, and a nearly expo- 
nential cut-off at a characteristic energy E c of a few hundred 
keV (see e.g. Gierlinski et al. 1997). Occasionally, the source 
switches to the High Soft State (HSS). The high-energy power- 
law is then much softer (V > 2.4) and the bolometric luminos- 
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Fig.l. RXTE/ASM light curve of Cygnus X-l. The time of 
the INTEGRAL observation is delimited by the vertical dashed 
lines 



ity is dominated by a thermal component peaking at a few keV. 
Finally, there are also Intermediate States (hereafter IMS) in 
which the source exhibits a relatively soft hard X-ray spectrum 
(F ~ 2.1-2.3) and a moderately strong soft thermal component 
(Belloni et al. 1996; Mendez & van der Klis 1997). The IMS of- 
ten, but not always, appears when the source is about to switch 
from one state to the other. When it is not associated with a 
state transition, it is interpreted as a 'failed state transition'. 
Until 1998, the source used to spend nearly 90 % of its time 
in the LHS. In the recent years however there have been more 
IMSs and soft states (see Zdziarski et al. 2002, Pottschmidt et 
al. 2003, Gleissner et al. 2004). 

Simultaneous radio/X-ray and high-energy observations of 
Cygnus X-l and other sources have shown that the X-ray LHS 
is correlated with a strong radio emission which is consistent 
with arising from a jet (Fender 2001). In contrast, during HSS 
episodes the source appears to be radio weak (Brocksopp et al. 
1999). The presence of a compact jet in the LHS was confirmed 
by Stirling et al. (2001) who presented evidence for an ex- 
tended and collimated radio structure on milliarsecond scales. 

State transitions are generally interpreted as being associ- 
ated with changes in the geometry of the accretion flow. In 
the HSS the geometrically thin optically thick disk (Shakura & 
Sunyaev 1973) extends down to the last stable orbit. The spec- 
trum is dominated by the thermal disc component and peaks 
at a few keV. The hard X-ray emission is then believed to be 
produced in a non-thermal corona above and below the disc 
(Gierliriski et al. 1999, hereafter G99). In the LHS the geo- 
metrically thin disc is truncated at a few hundred Schwarschild 
radii, the innermost part of the accretion flow forms a geomet- 
rically thick optically thin and hot disc (Shapiro et al. 1976; 
Narayan & Yi 1994) where high-energy radiation is produced 
trough thermal Comptonisation. During spectral transitions to 
the HSS the inner radius of the cold accretion disk decreases. 
This reduction of the inner disc radius is associated with either 
the cold disk penetrating inside the hot inner flow, or the later 



collapsing into an optically thick accretion disk with small ac- 
tive regions of hot plasma on top of it (Zdziarski et al. 2002). 
In both cases the enhanced soft photon flux from the disk tends 
to cool down the hot phase, leading to softer spectra. An al- 
ternative possibility would be that in both LHS and HSS, the 
geometry is that of a corona above a cold standard accretion 
disc (Bisnovatyi-Kogan & Blinnikov 1976; Haard & Maraschi 
1993). In the LHS the coronal plasma is essentially thermal 
and the cold disc faint because most of the power is trans- 
ported away and dissipated in the corona and a strong outflow 
(Beloborodov 1999; Malzac, Beloborodov & Poutanen 2001; 
Merloni & Fabian 2002; Ferreira et al. 2005), while in the HSS 
most of the power is dissipated into the cold disc and the corona 
is non-thermal. 

In the LHS, the radio flux is then positively correlated with 
the soft X-ray emission (3 - 25 keV, Corbel et al. 2000, 2003; 
Gallo, Fender, Pooley 2003). The X-ray emitting region and 
the jet seems to be intrinsically associated. This led to the now 
widely accepted idea that the corona (or hot thick disc) of the 
LHS constitutes the base of the jet (Fender et al. 1999; Merloni 
& Fabian 2002; Markoff, Nowak & Wilms 2005) 

Cygnus X-l represents a prime target for the INTEGRAL 
mission (Winkler et al. 2003) launched in 2002 October 17, 
whose instruments offer an unprecedented spectral coverage at 
high-energy, ranging from 3 keV to several MeV. Cygnus X-l 
was extensively observed (1 Ms) during the calibration phase 
of the mission (Bouchet et al. 2003, Pottschmidt et al. 2003, 
Bazzano et al. 2003). At that time, the source presented all the 
characteristics of the LHS. The source was later observed again 
during the galactic plane survey and open-time programme, 
and occasionally found in IMSs (Cadolle Bel et al., 2005, here- 
after CB05). In this paper we focus on the results of the first ob- 
servation of Cygnus X-l in the open time programme. This 300 
ks observation was performed on 2003 June 7-11 (rev 79/80) 
with a 5 x 5 dithering pattern 1 (the effective exposure time 
was 275 ks for JEM-X2, 292 ks for IBIS/ISGRI, and 296 ks 
for SPI). At this epoch, the RXTE All Sky Monitor count rate of 
Cygnus X- 1 was higher than in typical LHS by up to a factor of 
4, and the light curve showed strong X-ray activity characteris- 
tic of state (or failed state) transitions (see Fig.[Q. We also com- 
bine the INTEGRAL data with the results of coordinated radio 
observations (15 GHz) performed with the Ryle telescope. In 
Sec. |3 we present a spectral analysis of the JEM-X, IBIS/ISGRI 
and SPI spectra averaged over the whole duration of the obser- 
vation, in Sec.|3we study the strong broad band variability of 
the source during the observation. 

2. Average Spectrum 

2.1. Data processing 

We reduced the IBIS/ISGRI and JEM-X data with the standard 
analysis procedure of the Off-Line Scientific Analysis OSA 

4.2. released by the ISDC, whose algorithms are described 



1 INTEGRAL observations are made of a succession of exposures of 
about 30 minute duration with varied pointed directions to enable SPI 
image deconvolution. Such a 30 minute pointing is called a science 
window 
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Fig. 2. Best fits of the IBIS/ISGRI (red) and SPI (blue) data with a pure comptonization model (compps, left) and Comptonisation 
plus power-law (right). 
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Fig. 3. Joint JEM-X/SPI/ISGRI spectrum of Cygnus X-l av- 
eraged over revolutions 79 and 80. The data are fitted with 
the thermal/non-thermal hybrid Comptonisation model eqpair 
with power-law injection of relativistic electrons (see text 
and table |3. The lighter curves show the reflection com- 
ponent (solid), the disc thermal emission (dashed) and the 
Comptonised emission (dot-dash). The green, red and blue 
crosses show the JEM-X, IBIS/ISGRI and SPI data respectively 



in Goldwurm et al. (2003) and Westergaard et al. (2003) for 
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Fig. 4. Same as in Fig. except that the fitting model is now 
eqpair with mono-energetic injection of relativistic electrons 
with Lorentz factor y; n j =8.6 (see text and tabled 



IBIS and JEM-X respectively. A basic selection was performed 
to exclude those pointings too close to radiation belt entry or 
exit or spoilt by too large noise. The SPI data were prepro- 
cessed with OSA 4.2 using the standard energy calibration 
gain coefficients per orbit and excluding bad quality pointings 
which have anomalous exposure and dead time values, or with 
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Fig. 5. Same as in Fig. [3] except that the fitting model is now 
eqpair with a power-law injection relativistic electrons plus an 
additional warm Comptonisation soft component (comptt, see 
text and table|2jl. The comptt component is shown by the 3-dot- 
dashed curve 



a high final chi-squared during imaging, spiros 9.2 (Skinner & 
Connell 2003) was used to extract the spectra of Cygnus X- 
1, Cygnus X-3 and EXO 20390+375, with background model 
proportional to the saturating event count rates in the Ge de- 
tector. Concerning the instrumental response, version 15 of 
the IRF (Image Response Files) and version 2 of the RMF 
(Redistribution Matrix Files) were used. 

We produced JEM-X (3-25 keV), ISGRI (20-800 keV) and 
SPI (25-1000 keV) spectra averaged over revolution 79 and 80. 
Uncertainties of 3 % were added in quadrature to all three spec- 
tra to account for systematic errors. The resulting spectra were 
fitted using xspec v 11.3.1 

2.2. JEM-X spectrum 

The best power-law fit to the 3-20 keV JEM-X spectrum has 
a spectral index T = 2.16 which is intermediate between LHS 
and HSS spectral indices. The fit is however not statistically 
acceptable (x 2 /v = 611/131) with strong residuals indicating 
the presence of strong reflection and broad iron line features 
and a soft component (see below in Sect. 12. 41 and Fig. [3}. 

2.3. SPI and IBIS/ISGRI spectra 

The spectra show a highly significant evidence for a cut-off (or 
at least a break) around 100 keV (see Fig. |2jl. We first fit the 
ISGRI and SPI spectra independently with the Comptonisation 
model compps (Poutanen & Svensson 1996) The compps model 
provides a reasonable fit with reduced^ 2 values 0.82 and 1.12 



Table 1. Best-fit parameters of the SPI and IBIS spectra fit- 
ted with thermal-Comptonisation and thermal-Comptonisation 
plus power-law models. For each instrument the first line gives 
the results of the compps model alone (Comptonisation tamper- 
ature kT and Thomson depth r) while the second line gives the 
results of the fit with compps plus a power-law (photon index 
T). In the compps model the black body temperature of the soft 
seed photons was fixed at 0. 1 keV in all fits. 
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SPI+ISGRI 


mi 


i nn+0.14 
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SPI+ISGRI 


42!> 5 


2 64 +0 - 50 


9 9*7+0.66 
z - z '-0.77 


106/100 



for ISGRI and SPI respectively. The best-fit parameters ob- 
tained for both instruments (shown in table ^ are compatible 
within the error bars. 

However the residuals show an excess emission with re- 
spect to the Comptonisation model above the thermal cut-off 
in the SPI spectrum, which is not detected in the ISGRI spec- 
trum. This prompted us to reprocess the SPI data using differ- 
ent background models based on an empty field observation to 
model the uniformity map of the detection plan. We also tested 
using the instrument team processing software instead of OS A. 
This did not affect the presence of a high energy excess in the 
SPI spectrum. We concluded that the non detection of this ex- 
cess by ISGRI was attributable to its poorer sensitivity above 
200 keV. This is also the reason why the ISGRI fit appears sta- 
tistically better. Within the errors bars the highest energy IBIS 
and SPI data points are not too far and still compatible. 

To obtain better constraints on the parameters we then fit 
simultaneously the SPI and ISGRI data. The results are shown 
in Fig. |2 In all fits we allow for a free normalisation constant of 
the model for each instrument to correct for inter-calibration er- 
rors. The resulting difference in normalisation between the in- 
struments never exceeds 20 %. The simple compps model leads 
to best-fit parameters that are intermediate between those ob- 
tained for the individual ISGRI and SPI spectra, but the sta- 
tistical quality of the fit is not as good (x 2 /v = 1.4). As a 
further test for the presence of a high-energy excess, we re- 
fit both spectra with compps adding a power law component. 
When we allowed for this additional component the reduced 
X 1 decreased to 1.05. This improvement is highly significant, 
a F-test (Ax 2 = 36 for 2 additional parameters) shows that 
the probability that this improvement occurred by chance is 5 
xl0~ 7 . We conclude that either the high-energy cut-off is not 
due to thermal Comptonisation (e.g. associated to a cut-off in a 
non-thermal lepton distribution) or there is an additional non- 
thermal component at high-energy. We stress that the presence 
of the high energy excess is not in conflict with the ISGRI data 
although this instrument was not able to detect it. 
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Table 2. Best-fit parameters of the joint JEM-X, SPI 
and IBIS/ISGRI spectra with hybrid thermal/non-thermal 
Comptonisation models (eqpair, see Sect. 12.41 . The temper- 
ature of the inner disk (diskpn) was fixed to kT mlix =0.3 keV 
in all fits. The soft photon compactness is fixed at l s = 10. 
The absorbing column density is Nh = 5 x 10 21 and the in- 
clination angle 45 degrees. The fit parameter r p refers to the 
Thomson optical depth of the electrons associated with ions. 
The table also gives the resulting total optical depth tt, in- 
cluding electron-positron pairs whose production is calculated 
self-consistently, and the temperature kT e of the thermalised 
particles, computed according to energy balance. The extrapo- 
lated 0.1-1000 keV model flux, Fb i, and thermal disc compo- 
nent flux, are given in units of 10~ 8 ergs s _1 cirT 2 . The 
first column shows the results of the fit with a power law in- 
jection of non-thermal electrons with Lorentz factors ranging 
from 7min=1.3 to y m ax=1000. The data are then fitted for the 
index of the distribution of injected electrons T p . The second 
column shows the best fit parameters for a mono-energetic in- 
jection of electrons, the fit parameter y;„j represents the Lorentz 
factor of the injected particles. Finally, the third column gives 
the results for a power-law injection plus an additional warm 
Comptonisation component (comptt). Its best fit temperature 
AT comptt , Thomson depth r compt t and flux .Fcomptt are also shown. 
In this fit F p pegged to its minimum allowed boundary (2) 



Model 



T v 

T p or y inj 

£ (erg cms" 1 ) 
£i inc (keV) 
EW (eV) 
kT comptt (keV) 

^"comptt 
^comptt 

kT s (keV) 

tt 
* 2 M.o.f. 



EQPAIR pOW. 
1 1Q+U.U7 
i - i? -0.06 
1+0 

-0.04 
i t | +0.02 
^ -0.04 

0.67i°;» 
2063^ 3 
6 76 +0 - 36 

-0.33 

55 +33 

-"-34 



eqpair mono. 

0-51^4 
55 +0 01 
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8 41+ 62 
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u. / ^_ 003 
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7 02 +0 - 32 

' '"^-0.23 
90+38 
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0.79+° 



-0.1)2 
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9+0.21 
A -0 
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u ' u -0.06 
qq 1 +297 
yyl -176 

6 67 +0 25 



76.3!§ 

9 co+0.20 
—0.99 
9 7T+0.34 
-0.27 







1.52 


1.9 


3.0 


0.63 


7.5 


7.9 


5.9 


20.3 


50 


65 


1.36 


0.55 


1.01 


316/245 


244/245 


216/242 



2.4. Joint spectrum 

We now consider the joint JEM-X/ISGRI/SPI spectrum. The 
simple models considered in the previous section do not pro- 
vide an acceptable description of the broad band spectrum: 
strong reflection features and a soft-excess are required by 
the JEM-X data. CB05 present spectral fits of the same data 
with various thermal and non-thermal Comptonisation mod- 
els. In these models the soft excess in the JEM-X spectrum 
is accounted for by thermal emission of the accretion disc. 
However, although, these fits are statistically acceptable, they 
require high temperatures of the accretion disc (1.21+0.29 
keV). As mentioned in CB05 such a high temperature of the 
accretion disc is not realistic. Since the distance and scale of 
the accretion disc (i.e. mass of the black hole) are quite well 



known for this source, the intrinsic luminosity of the accretion 
disc is constrained: if the accretion disc in Cygnus X-l had an 
inner temperature of 1 keV, it would produce a flux at least one 
order of magnitude larger than what is observed. 

In order to get a more physically motivated representa- 
tion of the data, we investigate in more details thermal/non- 
thermal hybrid Comptonisation models and attempt to fit the 
data using different variants of the eqpair model (Coppi 1999; 
G99;Fronteraetal. 2001 (hereafterF01);Zdziarskiet al. 2002, 
2004) where the accretion disc temperature is arbitrarily fixed 
to a reasonable value. This model was shown to be successful 
in accounting for the high energy spectra of Cygnus X-l and 
other black holes candidates in different spectral states and over 
a very broad energy band ranging from soft X-rays to gamma- 
rays ( see e.g. McConnell et al. 2000, 2002). The wide use of 
this model in the literature will allow us to compare our obser- 
vation with previously published results. 

2.4.1. The hybrid thermal/non-thermal Comptonisation 
model (eqpair) 

A detailed description of the eqpair model can be found in G99. 
Its main ingredient is a spherical hot plasma cloud with contin- 
uous acceleration of electrons intended to model the emission 
of the hot disc/coronna. The high-energy electrons lose energy 
because of Compton, Coulomb, and bremsstrahlung processes 
and thus establish a steady-state distribution. At high energies, 
the distribution is non-thermal (power-law like), but at low en- 
ergies a Maxwellian distribution is established. The tempera- 
ture of the Maxwellian population, kT e is determined by bal- 
ance between Compton gains and losses, Coulomb heating by 
high energy electrons, bremsstrahlung losses, and direct heat- 
ing (e.g., Coulomb heating by energetic ions). The total num- 
ber of electrons (not including e + - e~ pairs, the production of 
which is also taken into account) is determined by the corre- 
sponding Thomson optical depth t p which is a free parameter. 
The cloud is illuminated by soft thermal photons emitted by an 
accretion disk. These photons serve as seed for Compton scat- 
tering by both thermal and nonthermal electrons. The system 
is characterised by the power L, supplied to its different com- 
ponents. We express each of them dimensionlessly as a com- 
pactness, lj = L,a"r/(%WeC 3 ). % is the characteristic dimension 
of the plasma, cr T is the Thomson cross-section. l s , l t h, /„,/,, and 
h- Uh+lnth correspond to the power in soft disk photons enter- 
ing the plasma, thermal electron heating, electron acceleration 
and the total power supplied to the plasma. We follow G99 and 
fixed l s = 10. 

The disc spectrum incident on the plasma is modelled 
as coming from a pseudo-Newtonian accretion disk extend- 
ing from R out = 10 3 R g down to the minimum stable orbit, 
Ri n — 6R g . Its spectral shape is then characterised by the max- 
imum colour temperature of the disk, kT max . Previous obser- 
vations of Cygnus X-l with X-rays telescopes indicate tem- 
peratures ranging from 0.1 keV in the LHS to up to 0.6 keV 
in the HSS (G99;F01). For such temperatures the peak of the 
thermal disk emission is below the energy range covered by 
INTEGRAL, and the disk temperature cannot be constrained by 
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our observations, we therefore fixed kT ma x = 0.3 keV. The cov- 
ering factor of the corona is unity. 

The spectrum from both reflection and the Fe K a is calcu- 
lated, taking into account relativistic smearing with the emis- 
sivity dependence oc R 2 . The reflecting material is allowed 
to be ionised, with the degree of ionisation characterised by 
the ionisation parameter We follow G99 and fixed the col- 
umn density of absorbing material along the line of sight to 
Nh — 5 x 10 21 crrr 2 and the inclination angle of the system at 
45 degrees. 

2.4.2. Fit results 

We first performed a fit similar to those of G99 and F01 with the 
non-thermal electrons injected with a power-law distribution of 
Lorentz factors ranging from y m \ n = 1.3 to y max = 1000. The 
upper and lower limits y m ; n and y max where kept fixed while 
fitting for the power-law index F p . This results in an accept- 
able fit with a reduced^- 2 of 1.29. The unfolded broadband 
spectrum and residuals are shown in Fig. [3] The best-fit pa- 
rameters are presented in Table 13 1^1 l s is about unity i.e. in- 
termediate between what is generally found in the LHS (4-10, 
see Ibragimov et al. 2005) and the HSS (< 0.4), the heating 
of the plasma is dominated by the non-thermal acceleration 
{Imhlh = 1)- Overall the parameters are similar to those ob- 
tained by G99 for the IMS of 1996 May 23. However an inspec- 
tion of the residuals (Fig. [5} shows that the high energy part of 
the ISGRI/SPI spectrum is not well represented by this model. 
The model does not provide a good description of the shape 
of the high-energy cut-off and overestimates the measured flux 
above 200 keV. As an alternative we attempted to fit the spec- 
trum assuming that the electrons are injected at a single Lorentz 
factor yi n j instead of a power-law distribution. Such a mono- 
energetic injection is not expected in the case of shock accel- 
eration but could be achieved in reconnection events that are 
expected to power the corona. The resulting fit is displayed on 
Fig.0]and the unabsorbed best fit model spectrum is shown on 
Fig-El The best-fit parameters are shown on Tabled Note that 
the fraction of non-thermal power is now only - 0.51. 

This model gives a good description of the INTEGRAL spec- 
trum. The best-fit value for the maximum Lorentz factor of the 
electrons y; n j =8.6 implies that the once scattered soft pho- 
tons reach a maximum energy of about 4y 2 a -kT m!a ~ 90 keV. In 
other words, in this model the high-energy cut-off observed in 
the spectrum is non-thermal and corresponds to the upper end 
of the non-thermal electron distribution. 

So far we have shown that the soft excess in the JEM-X 
data does not require an unphysical high accretion disc temper- 
ature and can be accounted for by an hybrid particle distribu- 
tion. In this context, our data favour mono-energetic injection 
over power-law acceleration. However, this is not the only pos- 
sible interpretation. Indeed, such a soft excess is not uncom- 
mon in IMSs and the softest LHSs. It is often interpreted as a 
hot spot on the accretion disc or alternatively as a component 
due Comptonisation in the warm upper layers of the disc. It can 
be accounted for by adding a second Comptonisation compo- 
nent with a temperature of a few keV (Di Salvo et al 2001 ; F01 ; 



Zicky, Done & Smith 2001). Alternatively, Markoff, Nowak & 
Wilms (2005), suggest that this soft component originates from 
optically thin synchrotron emission in the jet. 

To check the effect of such an additional soft component 
on the derived parameters of the eqpair model, we fit the spec- 
trum with a model consisting in the power-law injection eq- 
pair model, including the same 0.3 keV thermal component, 
relativistically smeared iron line and ionised reflection compo- 
nent as in our previous models (see Sect. 12.4. 1> plus an addi- 
tional Comptonisation component (comptt model in XSPEC) 
in a warm medium. 

This led to an excellent fit of the INTEGRAL data, with 
X 2 /y = 0.89, that is shown in Fig. |5] and Table [2] The addi- 
tional Comptonisation component has a best-fit temperature of 
2.5 keV and a Thomson optical depth of 2.7. The fit is much 
better than what obtained with the powerlaw injection eqpair 
model. The probability that the improvement obtained when 
the soft comptt component is added to the powerlaw injection 
eqpair model occured by chance is 7x 10~ 20 according to F- 
test (Ay 2 = 100 for 3 additional parameters). This fit is statis- 
tically comparable to that obtained with mono-energetic injec- 
tion. It demonstrates that mono-energetic injection is far from 
inevitable. 

We also note an interesting difference between power-law 
and mono-energetic injection models: the presence of a broad 
e + —e~ pair annihilation line around 511 ke V that is clearly visi- 
ble in Fig.[3]and |5](power-law injection) but absent from Fig.0] 
(mono-energetic injection). In the power-law model some elec- 
trons are injected at very high Lorentz factors (< 1000) these 
electrons efficiently up-scatter soft-photons at energies above 
the pair production threshold, leading to a pair cascade and 
subsequent annihilation, whereas in the mono-energetic case 
the pair content of the plasma is negligible due to the much 
lower energy of the electrons (< 8.41). Unfortunately the poor 
statistics of the data does not allow us to discriminate between 
these models on the basis of the annihilation line. 

In the eqpair+comptt spectral model the fraction of non- 
thermal power (Znth/'h) is reduced to about 20 % (to be com- 
pared to 100 % and 50 % obtained respectively for the simple 
power-law and mono-energetic injection models). This shows 
that the measured non-thermal fraction is sensitive to the as- 
sumptions of the spectral model. In the latter fit the hot plasma 
is essentially thermal, in agreement with the conclusions of 
CB05. Indeed, the equilibrium temperature and optical depth 
we obtain are similar to what obtained by these authors fitting 
the same data with a simple thermal Comptonisation model 
plus a multicolour disc spectrum. We also note that in the eq- 
pair+comptt fit, the hard to soft compactness ratio is larger 
(lh/k — 3 instead of ^ 1 in the previous fits). Overall the pa- 
rameters are closer to what is usually obtained in the LHS. 

In short, our spectral analysis does not favour the simple 
power-law injection eqpair model. Rather, models including 
an additional soft component (and a weak non-thermal frac- 
tion), or alternatively, models with an important thermal frac- 
tion but mono-energetic injection are preferred. There are how- 
ever some caveats in our spectral analysis that might affect any 
conclusion drawn from our fits: First, there are still open cali- 
bration issues, in particular regarding the X-ray monitor JEM- 
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X specially below 5 keV. It is in this energy range that the soft 
component shows up and it is very important to constrain the 
parameters of the hybrid model. Second, our model assumes a 
unique emitting zone although it is likely that in an IMS the 
situation is more complex. Indeed, the thermal hot flow of the 
LHS probably coexists with the non-thermal corona of the HSS 
in distinct regions of the accretion flow. If various emitting re- 
gions, with very different spectra and physical parameters, con- 
tribute to the observed spectrum, then the derived best-fit pa- 
rameters have a poor physical significance. Similarly, and per- 
haps more importantly, the parameters determining the shape 
of the spectrum vary not only within the accretion flow but also 
with time. As we will show below, the source exhibited a strong 
spectral variability during the observation with the appearance 
of spectra ranging from quasi-LHS to quasi-HSS. Under such 
circumstances, the precise value of the best-fit parameters of 
the average spectrum may be physically irrelevant. 

3. Spectral Variability 

3.1. Setup 

In order to study the spectral variability of the source during the 
observation, we produced light curves in 16 energy bands rang- 
ing from 3 to 200 keV 2 with a time resolution of the duration 
of a science window (i.e. ~ 30 min). Above 200 keV, the vari- 
ability is dominated by statistical noise. The count rate in each 
band was then renormalised so that its time average matches 
the energy flux calculated from the best-fit model of the joint 
average JEM-X/ISGRI/SPI spectrum shown in Fig. [5] Namely, 
for each energy band we compute the quantity 



C(t) - 
F(f) = -=-F, 
C 



(1) 



where C(f) is the mean count rate during pointing t, C is the 
count rate averaged over the whole observation, F is the ob- 
servation average energy flux in this band given by the best 
fit model. We use F(t) as a proxy for the instantaneous en- 
ergy flux, therefore neglecting the effects of the spectral vari- 
ations on the instrumental response. This simple deconvolu- 
tion method provides us with a convenient approximate of the 
broadband energy spectrum for each pointing that will be use- 
ful for a physical interpretation of the variability. This method 
is much more convenient than fitting the spectra for each sci- 
ence window, in particular if one considers the large number 
of parameters required to fit the data and the poor statistic in 
the short exposure spectra. On the other hand our method for 
estimating the spectra from the light curves also enables us to 
improve the photon statistics by combining the IBIS/ISGRI and 
SPI instruments in the energy range where they overlap. The 
time dependent flux F(t) in an overlapping band is estimated 
as follows: 



F(t) = 



of Fi(0 + a\ F s (t) 



erf + of 



(2) 
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Fig. 6. Light curves in different JEM-X an ISGRI/SPI bands 
with increasing photon energy from bottom to top, as labelled 
in keV. Each light curve was re-scaled to the average flux and 
then incremented according to its energy range for clarity. 



2 namely: 3^1, 4-5, 5-6, 6-7, 7-9, 9-11, 11-13, 13-15, 15-20, 20- 
30, 30-40, 40-50, 50-80, 80-100, 100-140 and 140-200 keV 



Where F\ and Fs are the measured ISGRI and SPI fluxes in 
that band; and cr s are their time averaged uncertainties. This 
combination minimises the average uncertainty on F. 

The resulting light curves are shown in Fig. [6] The time 
averaged 3-200 keV model flux is F3-200 = 2.87 x 10~ 8 ergs 
cm" 2 s . The energy flux has a rms amplitude of 16 %, and the 
ratio of the maximum to the minimum luminosity is 2.6. 
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Fig. 10. PCA parameters a\ (left panel) and ai (right panel) as a function of the radio flux (diamonds). In both panels, the best 
linear fits are shown by the solid lines. The crosses indicate the time average over each of the five periods of nearly continuous 
radio coverage (see Fig |9j- While there is no convincing correlation between the radio flux and a\, the radio flux is correlated to 
<X2 at highly significant level. 



3.2. Principal Component Analysis (PCA) 

The light curves shown in Fig.|6]exhibit a complex and strong 
broad band variability of the spectra as well as the overall flux. 
We use a principal component analysis (PCA) to seek for vari- 
ability patterns in our sample. PCA is a powerful tool for mul- 
tivariate data analysis used for a broad range of applications in 
natural as well as social sciences (see e.g. Kendall, 1980). It 
has also been used in astronomy. For instance, Francis & Wills 
(1999) provide a brief introduction to PCA as applied to quasar 
spectra. Previous application of PCA to spectral variability in- 
clude Mittaz, Penston & Snijders (1990) who discussed the UV 
variability of NGC 415 1 and Vaughan & Fabian (2004) for the 
X-ray variability of MCG6-30-15. 

The main use of PCA is to reduce the dimensionality of 
a data set while retaining as much information as possible. 
It transforms a number of (possibly) correlated variables into 
a (smaller) number of uncorrected variables called principal 
components. These principal components may define patterns 
or correlations in the data set that can often be interpreted more 
easily than the original (large) data sets. 

In our specific application we have p — 79 spectra mea- 
sured at times t\, to, ..„t p and bined into n — 16 bins cor- 
responding to energies E\, E%,...,E n . Each spectrum can be 
thought of as a point in a n-dimensional space, so that the co- 
ordinate of the point j along the £, axis is given by the cor- 
responding energy flux F(tj, £,). The shape and extension of 
the cluster of points formed by the whole set of spectra charac- 
terises the source variability. The idea behind PCA is to deter- 
mine a new coordinate system in which the description of this 
cluster will be simpler. 

For this purpose, the data can be viewed as the p x n matrix 
which coefficients are given by the energy fluxes F(tj, Ef). Then 
one can compute the n x n covariance matrix of the data. In 



practice, PCA consists in diagonalizing this covariance matrix 
to obtain its eigenvalues and eigenvectors. This procedure gives 
the coordinates of each eigen vector i.e. the Ck(Ei) coefficients 
denoting the coordinate of the fc-th eigenvector along £, axis. 

These eigenvectors define a new coordinate system, in this 
n-dimensional parameter space, which best describes the vari- 
ance in the data. The first principal component, or PC 1 (the 
eigenvector with the highest eigenvalue), marks the direction 
through the parameter space with the largest variance . The next 
Principal Component (PC 2) marks the direction with the sec- 
ond largest amount of variance, etc... 

Let aii(t j) be the new coordinate of the spectrum j along the 
k-th eigen vector, then the relation between the two coordinate 
systems can be written as follows: 

n 

F{t lt Ed = F{E,) + a k (tj)C k (Ed, (3) 

k=l 

where F(Ef) is the time averaged flux at energy £, (the time 
averaged spectrum is used as the origin for the eigenvector co- 
ordinate system). 

Eq.|3]amounts to a linear decomposition into n independent 
components of the variability (the eigenvectors C\ , C2,...,C„). 
The normalisation coefficients of each PCA component (re- 
spectively a\, ar2,..., a„) vary in time. Their fluctuations ac- 
count for the sample variance. On the other hand the eigenvec- 
tors Ck are constant, they define the variability mode of each 
PCA component. Since the eigenvectors describe fluctuations 
around the time-average flux, both the and Cu coefficients 
can take negative as well as positive values. If the Ck coeffi- 
cients are positive at all energies, the variability mode associ- 
ated to component k can be understood as due to an additive 
spectral component with a fixed shape and a variable normal- 
isation. When the sign of Q depends on energy, this corre- 
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Fig. 7. The 2 first principal components of variability. The up- 
per panels illustrate the effects of the each component on the 
shape and normalisation of the spectrum: time average spec- 
trum (light line) and spectra obtained for the maximum (thick 
solid line) and minimum (dashed line) observed values of the 
normalisation parameter. The middle panels show the ratio of 
spectra obtained of the maximum and minimum spectra to the 
average one. The bottom panels show the contribution of each 
component to the total variance as a function of energy. 



sponds to more complex spectral variability modes (e.g. pivot- 
ing). 

As mentioned above, the PCA components are ordered ac- 
cording to the amount of sample variance they account for (i.e. 
the observed fluctuations of a\ cause more variance than those 
of ff2 which produce more variance than 0-3 etc.). The first 
few Principal Components (those representing most of the vari- 
ance in the data) should reveal the shape of the relevant spec- 
tral components or variability modes. The weaker Principal 
Components might be expected to be dominated by the statis- 
tical and systematic noise in the spectra. To summarise, PCA 
finds the decomposition that maximises the variability due to 
lower order components, so that most of the variability can be 
described using a small number of components. 




"5 1.00 
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Fig. 8. Observed variance spectrum (squares). The curves show 
the contribution of PC 1 (3 dots-dash), PC 2 (dashes), PC 
3 (dot-dash), estimated statistical noise (dots) and their sum 
(solid). The bottom panel shows ratio of the contribution of 
PC 1+PC 2+PC 3+noise to the total observed variance (values 
exceeding unity might be due to overestimated statistical noise) 

3.3. Results of the PCA analysis 

The results of our PCA analysis of the spectral variability of 
Cygnus X-l are illustrated in Fig.0 which shows how the 2 
first principal components affect the flux and spectrum and their 
respective contribution to the total observed variance as a func- 
tion of energy. As can be seen on this figure, the first principal 
component (PC 1) consists in a variability mode dominated by 
variations in the luminosity (normalisation) with little change 
in the spectral shape. For this reason, in the following, we will 
refer to PC 1 as the 'flaring mode'. This component accounts 
for 68 % of the sample variance. Let us now consider the small 
spectral fluctuations induced by PC 1 . An increase in luminos- 
ity is associated with spectral hardening in the 3-10 keV band, 
and a slight softening at higher energies which is particularly 
marked in the 10-30 keV band. This spectral variability sug- 
gests a variable power-law spectrum with fixed spectral index 
moving on top of a constant thermal emission disc emission 
plus reflection component. Despite these small spectral varia- 
tions, PCI correlates very well with the high-energy flux. A 
least square fit shows that the 3-200 keV flux relates to a\ 
through: 



^3-200 = ^3-200(1 + 7.64a!), 



(4) 



with a linear correlation coefficient of 0.98. So that a\ can be 
viewed as a tracer of the hard X-ray luminosity of the source. 

As shown in Fig0 the second PCA component (PC 2) can 
be described roughly as a pivoting of the spectrum around 10 
keV. The two spectra obtained for the minimum and maximum 
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Fig. 9. Evolution of the parameters associated to PC 1 (top), 
PC 2 (middle) and radio light curve (bottom) during the obser- 
vation 



values of the a<x parameter controlling the strength of PC 2 are 
reminiscent of the canonical LHS and HSS spectra. This com- 
ponent is responsible for 27 % of the sample variance, and will 
be referred as the 'pivoting mode', can thus be seen as a 
tracer of the hardness of the high-energy spectrum. 

The third PCA component (PC 3) accounts for only 2 % 
of the sample variance. It is dominated by fluctuations of the 
relative normalisation of the spectrum in the JEM-X band and 
ISGRI/SPI energy range. This component is most likely an in- 
strumental artefact due to known calibration issues related to 
the dithering observation mode. Indeed, since the pointed di- 
rection changes between successive science windows and as 
the effects of vignetting and background non-uniformity are not 
perfectly corrected in the present release of the data analysis 
software, this results in spurious variability that affects mainly 
the relative flux normalisation between instruments. This effect 
is particularly strong in JEM-X for which systematic errors as 
large as 30 % are expected for a source 5 degrees off axis. PCA 
enables us to disentangle such instrumental effects from the in- 
trinsic variability. Because they are independent of the intrinsic 
source variations and as long as they are not dominant, they are 
filtered and isolated into higher order PCA components. 

Similarly variability due to statistical noise is filtered out 
into higher order components. From the statistical uncertain- 
ties on the measured count rates we expect 3 % of the observed 
variance to be due to statistical noise. As the first 3 PCA com- 
ponents already account for 97 % of the sample variance, the 
higher order components are most likely due to noise. Fig. [8] 
shows the contribution of PC 1, PC 2, PC 3 and statistical 
noise to the observed variance spectrum (i.e. the measured vari- 
ance as a function of photon energy). These 4 components are 
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Fig. 11. Correlation between radio flux density and (100-140 
keV)/(4-5 keV) hardness 



enough to account for the observed variance within a few per- 
cents at all energies. 

Since PC 3 is likely to be an instrumental artefact and 
higher order components are probably noise, we conclude that 
the intrinsic source variability is largely dominated by PC 1 
(flaring mode) and PC 2 (pivoting mode). Fig.|5]shows the time 
evolution of the PCA parameters a\ and a%. a\, which traces 
the changes in bolometric luminosity at nearly constant spec- 
tra shows important variability on time scales of order of a few 
hours or less, but no clear systematic trend during the 4 days 
of observation. In contrast, az, which roughly traces the hard- 
ness of the spectrum, seems to vary on longer time scales: it 
jumps during the first 2 days then decreases in the second part 
of the observation. This suggests that the physical mechanisms 
responsible for PC 1 and PC 2 are distinct (which is also ex- 
pected from the fact that, by construction, PC 1 and PC 2 are 
linearly independent) and apparently acting on different time 
scales. 



3.4. Radio/high-energy correlation 

In order to study the possible correlations between the radio 
and hard X-ray emission, we selected the science windows for 
which we had simultaneous radio pointings. When there were 
several radio pointings falling within a given science window 
we use the mean radio flux. The resulting light curve is shown 
on the bottom panel of Fig. [9] A comparison with the time evo- 
lution of ct\ and aq. indicates that the radio flux tends to fol- 
low the evolution of the pivoting mode. In other words, the 
radio flux tends to be stronger when the hard X-ray spectrum 
is harder. This is clearly seen on Fig. EH which shows that the 
radio emission is strongly correlated to a<i. The correlation is 
highly significant. The Spearman rank test correlation coeffi- 
cient is 0.78 corresponding to a probability that the correlation 
is by chance of 2 x 10 7 . On the other hand, there is no hint 
of a correlation with the flaring mode as can be seen in the left 
panel of Fig.fTOI 
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Fig. 12. Flux (ergs crrT 2 s -1 ) in different energy bands (as indi- 
cated in keV), versus 15 GHz radio flux density 

It is also worth noting that the correlation between the piv- 
oting mode and the radio flux is valid only when we consider 
the data over more than one day, we find no evidence for a cor- 
relation on shorter time-scales. Therefore, our PCA analysis 
shows that on time scales of days, the radio jet activity is cor- 
related with hardness of the high-energy spectrum rather than 
hard X-ray luminosity. 

Obviously this can be seen directly from the light curves, 
although less comprehensively. For instance Fig. II II shows a 
strong correlation between the (100-140 keV)/(4-5 keV) hard- 
ness and the radio flux. Fig.ll2lshows that the radio flux tends 
to be anti-correlated with the X-ray flux (3-7 keV) keV) and 
correlated with the soft gamma-rays (>15 keV). This depen- 
dence of the INTEGRAL/iadio flux correlation confirms both 
the presence of a pivot point located around 10 keV and the 
correlation of the radio luminosity with the pivoting of the 
spectrum. However, without PCA it would have been very dif- 
ficult to demonstrate the presence of the 2 independent vari- 
ability modes. PCA also brings a much clearer information on 
the spectral variability due to each mode, and, on top of that, 
enables one to filter out systematic and statistical errors. 

4. Discussion 

We have shown that, during our IMS observation, the variabil- 
ity of Cygnus X-l can be described by two independent vari- 
ability modes: 

- On time-scales of a few hours or less there are important 
changes in luminosity with little spectral variations (flaring 
mode). 

- On longer time scales there is a spectral evolution with the 
the spectrum pivoting around 10 keV. 

We further showed that while there is no hint for a correla- 
tion between the radio flux and the flaring mode, the radio is 



strongly correlated with the pivoting of the spectrum, in the 
sense that the radio flux is stronger when the hard X-ray spec- 
trum is harder. This result strongly differs from what is usually 
reported in the LHS. Indeed, the radio flux is then positively 
correlated with the soft X-ray emission (3 - 25 keV, Corbel et 
al. 2000, 2003; Gallo, Fender, Pooley 2003). 

However, our results are not in conflict with these obser- 
vations since Cygnus X-l was not in a typical LHS but in an 
IMS, or rather was switching between different IMSs. Indeed, 
compilations of LHS and HSS spectra suggest that the spec- 
tral transition between LHS and HSS occurs through a pivoting 
around 10 keV (see e.g. Fig. 9 of McConnell et al 2002). The 
evolution of the parameter shown in Fig. [9] indicates that 
the source, initially in a 'soft' IMS, switched to a harder state 
during the first 2 days of observation and then transited back 
toward a 'soft' state. 

Actually, the transition from LHS to HSS is known to be 
associated with a quenching of the radio emission (Corbel et 
al. 2000; Gallo, Fender & Pooley 2003). As the transition to 
the HSS also corresponds to a strong softening of the spec- 
trum, this is consistent with the correlation between hardness 
and radio flux: when, during the observation, the source gets 
closer to the HSS the spectrum softens and simultaneously the 
radio flux decreases. We also note that a recent analysis of Ryle 
and RXTE data of Cygnus X-l (Gleissner et al., 2004) inter- 
estingly shows the same correlation tendencies during failed 
state transitions (Ryle/ RXTE/PCA: moderate anti-correlation, 
Ryle/ RXTE/HEXTE: correlation) as reported here, albeit on 
timescales from weeks to years. 

It is interesting to speculate on the cause of the two vari- 
ability modes. We tried to reproduce such variability modes 
by varying the parameters of the hybrid thermal/non-thermal 
Comptonisation models considered in Sect. 12.41 We used the 
best-fit model obtained with eqpair and mono-energetic parti- 
cle injection. We choose this model because it gives a better 
fit of the average INTEGRAL spectrum than the power-law in- 
jection model, still avoiding the complication of an additional 
component that cannot be computed self consistently. Actually, 
similar results are obtained for all of the 3 models we used to 
fit the broad band spectrum. 

As shown in the left panel of Fig.^]it is possible to pro- 
duce variations in luminosity by a factor comparable to what 
is observed and little spectral changes in the INTEGRAL band 
by varying the coronal compactness Zj, by a factor of 2. In this 
context the flaring mode would correspond to variations of the 
dissipation rate in the corona possibly due to magnetic recon- 
nection. This variability mode seems to be a characteristic of 
the HSS (Zdziarski et al. 2002). As we show here, it also pro- 
vides a major contribution to the variability of the IMS. 

Regarding the pivoting mode, it can be produced by 
changes in the flux of soft cooling photons at constant dissi- 
pation in the hot phase. We performed simulations assuming 
that the accretion disc radiates like a blackbody i.e. its flux 
Fdisc k Z s oc r^ ax and constant //,. For an increase of the disc 
temperature by a factor of 1.7, the disc luminosity grows by 
a factor of 8. As in this model, the disc flux also corresponds 
to the soft cooling photon input in the corona and the heating 
(oc l h ) is kept constant, this leads to a steepening of the spec- 
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Fig. 13. Left panel: effect of varying //, by a factor of 2 on the eqpair model with monoenergetic injection (see Sect. l2.4> . Solid 
curve: unabsorbed best-fit model (4 = 8.5); Dotted curve If, — 5.7; Dashed curve: //, = 11.9. Right panel: effect of varying the soft 
photons flux by a factor of 8. Solid curve: unabsorbed best-fit model (7di sc = 0.3 keV; lh/k = 0.85). Dotted curve: Tdi sc = 0.357 
keV and l h /l s = 0.42. Dashed curve: r d i sc = 0.212 keV and k/k = 3.4. 



trum with a pivot around 10 keV of similar amplitude as in PC 
2 (see Fig.Hl). For the 1996 HSS, G99 found a ratio l h /l s ~ 0.3 
while in the LHS lf,/l s ranges between 3.5 to 15 (Ibragimov et 
al. 2005). The range of Ihlh (0.4-3.4) required to reproduce 
the observed amplitude of the pivoting mode matches almost 
exactly the intermediate range between the HSS and the lower 
limit of the LHS. The source initially in a (quasi) HSS evolved 
toward the LHS but as soon as it was reached, it went back 
toward the HSS. 

Since, in the INTEGRAL band, the constraints on disc ther- 
mal emission are loose we did not attempt to model the data 
with a varying inner disc radius which is, moreover, difficult 
to disentangle from fluctuations of the disc temperature. In the 
fitted models as well as the models shown in Fig.^] the inner 
disc radius is fixed at 6 /?<,. Nonetheless our result would also be 
consistent with the disc moving inward and outward of the hot 
phase during the state transitions. Indeed, when the inner disc 
radius is approaching the black hole, its maximum temperature 
and luminosity increases 3 leading to a more efficient cooling 
of the hot flow/corona. The anti-correlation between radio flux 
and disc luminosity would be due to the jet expanding when 
the cold accretion disc recedes and then shrinking in the sec- 
ond phase of the observation when the disc moves back inward. 
It is interesting to note that the change in disc flux required (a 
factor of ~ 8) to explain the spectral evolution is comparable 
to the amplitude of the variations of the radio flux (a factor of 
~ 6). This suggests a direct relation between the disc flux and 

3 Unless the mass accretion rate is reduced by a larger amount, 
which seems very unlikely, the evidence being rather that the accretion 
rate is often (but not always) larger in the soft than in the LHS. 



jet power. The overall change in bolometric luminosity occur- 
ing during the PC2 transition estimated from the fiducial 'hard' 
and 'soft' state models shown on the left panel of Fig. is 
about a factor of 2. Because of the relatively short time scale 
(~ a day) on which the variation in luminosity occurs, it is un- 
likely to be driven by changes in the mass accretion rate. Most 
probably, it is due to a change in the radiative efficiency of the 
flow. The accretion flow could be less efficient in the LHS, be- 
cause about half of the accretion power is either swallowed by 
the black hole or pumped into the jet, while, in the HSS, the 
cold disc is expected to be radiatively efficient. 

The contribution of the jet power to the total energetic out- 
put of black hole candidates in the LHS is a matter of debate. 
It has been argued that the jet could be strong and even dom- 
inant (Fender et al. 2003; Malzac, Merloni & Fabian 2004). 
The jump in luminosity by a factor of 2 that we infer during 
our mini state transition sets an upper limit to the jet power, 
namely, the jet power at a luminosity just below the transition 
is at most comparable to the X-ray luminosity. This uppper 
limit is in aggreement with a recent study of the jet interaction 
with the surrounding interstellar medium by Gallo et al. (2005). 
According to these authors, the jet power in Cygnus X-l repre- 
sents between 6 to 100 % of the X-ray luminosity at the peak 
of the LHS. In our model, the upper limit is reached if there is 
no advection into the black hole. In this case the mini-transition 
pivoting mode would correspond to a redistribution of the ac- 
cretion power between the jet and the cold accretion disc. In 
this case all LHS sources would be jet dominated, since, ac- 
cording to the jet and X-ray power scaling laws of Fender et al. 
(2003), the jet share of the energy budget is increasingly larger 
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at lower luminosities. However, the amplitude of the luminos- 
ity changes we infer during the mini-transition (i.e. the factor of 
2) is model dependent and also depends on the strength of the 
thermal disc contribution to the time-averaged spectrum, which 
is poorly constrained with the INTEGRAL data. Similar obser- 
vations with and even broader spectral coverage (including soft 
X-rays) are required to consolidate this result. 

Moreover, although our data point toward an inefficient ac- 
cretion flow in the LHS, it does not tell us about the cause of 
this inefficiency: pure advection as well as jet dominated accre- 
tion flows are both viable possibilities. We note that a similar 
conclusion was reached by Chaty et al. 2003, on the basis of 
the analysis of the multi-wavelength spectrum of the transient 
black hole binary XTE Jl 1 18+480. 

The evolution of the hard X-ray corona luminosity dur- 
ing our IMS observation is very puzzling. Indeed, if, as com- 
monly believed for the LHS, the corona constitutes the base of 
the jet, it is difficult to conceive that changes in the jet power 
and/or extension is not associated to changes in the energet- 
ics of the corona. The apparent lack of response of the radio 
jet to the (short time-scale) X-ray fluctuations could well be 
due to the time delays required to propagate the information 
from the corona to the distant radio emitting region, which 
moreover have very different sizes. However, one would still 
expect the corona/hot accretion flow to track the longer time- 
scale evolution of the radio jet and also respond to changes in 
the disc power and/or distance of the truncation radius. Instead, 
we infer dramatic changes in the jet and disc power that are 
anti-correlated with each other, but completely unrelated to the 
fluctuations of the coronal power. On the one hand, the spectral 
pivoting described by PC 2 is understood in terms of changes 
in the disc luminosity at constant coronal luminosity. And, on 
the other hand, the rapid fluctuations of the coronal power that 
we do observe through PC 1 are apparently not associated to 
fluctuations of the cold disc or jet emission. The nature of the 
instabilities responsible for a coronal activity that is so uncou- 
pled to the jet and disc emission remains to be clarified. In any 
case, these results indicate that, in the IMS, the corona does not 
play the same role as in the LHS, and the whole disc-corona-jet 
interaction seem to work differently. 



5. Conclusion 

The INTEGRAL IMS spectrum of Cygnus X-l shows a high- 
energy cut-off or break around 100 keV. The shape of this 
cut-off differs from pure thermal Comptonisation, suggesting 
the presence of a non-thermal component at higher energies. 
The average broad band spectrum is well fitted with hybrid 
thermal/non-thermal Comptonisation models, although some 
important parameters such as the fraction of non-thermal power 
are not well constrained because they depend on the assump- 
tions of the model (mono-energetic versus power-law injection, 
presence or absence of an additional soft component). Models 
with mono-energetic injection, or models with an additional 
soft component seem to be favoured over standard power-law 
acceleration models. 



During our observation, the source presented a strong flux 
spectral variability occurring through 2 independent variability 
modes: 

i) changes in the dissipation rate in the corona, due to local 
instabilities or flares, producing a variability of the hard X- 
ray luminosity on time-scales of hours and no strong spec- 
tral alterations. Strikingly, this coronal activity seems to be 
unrelated to the evolution of the jet and cold disc luminos- 
ity. 

ii) a slower 4-day evolution starting from a spectrum close to 
the canonical HSS toward an almost LHS and back. This 
spectral evolution was characterized by a pivoting of the 
spectrum around 10 keV. It was correlated with the radio 
emission which was stronger when the hard X-ray spectrum 
was harder. It is interpreted in terms of a variable soft cool- 
ing photon flux in the corona associated with changes in the 
thermal disc luminosity and radio-jet power. This interpre- 
tation suggests a jump in bolometric luminosity of about a 
factor of 2 during the transition from LHS to HSS, which 
might indicates that the LHS accretion flow is radiatively 
inefficient, half of the accretion power being possibly ad- 
verted into the black hole and/or the radio jet. 

In the IMS, the jet power appears to be anti-correlated with the 
cold accretion disc luminosity, while the coronal power fluctu- 
ates independently. Apparently, the coupling between accretion 
and ejection processes differs from that of the LHS where the 
radio jet and the X-ray corona appear intrinsically linked. 
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